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(Abstract) The Fast Fourier transform (FFT) exist frequency domain spectral leakage, picket fence effect, as well as non- 
integer wave phenomenon. In this paper, a time-frequency filter is designed, which can detect the frequency, amplitude and phase 
of any order harmonics and interharmonics in signal by means of time domain convolution. Negative half-axis frequency is 
restrained by Hilbert transform. The theory analysis are carried to this method and the calculate formula are concluded. 
Experiment simulation results show that time -frequency filtering convolution flunction can be designed and realized neatly and 
be real-time implemented conveniently in engineering application; the influences of fundamental frequency fluctuation on 
harmonic analysis are restrained by using the approach presented in this paper;the calculating frequencies with many order 
harmonics and interharmonics are nicety, the amplitudes and the initial phases are detected in high-accurately. 
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1. INTRODUCTION 

High-precision analysis of harmonic power measurement, 
harmonic power flow calculation, network testing equipment, 
power system harmonics compensation and suppression is of 
great significance!!]. Because non-synchronous sampling and 
data truncation, the use of fast Fourier transform (FFT) 
algorithm to generate harmonic analysis and fence effect of 
spectrum leakage, the accuracy of harmonic analysis[2-3]. 

In order to eliminate these phenomena, scholars Design a 
lot of harmonic detection methods based on Fourier 
transform of many kinds window, such as rectangular [4], 
Hanning [5], Hamming [6], Blackman [7], Blackman-Harris 
[8], Kaiser window[9] and various improvements [10-11] and 
other windowed. Fast Fourier Transform(FFT) can reduce the 
encounter alone and fence effect of spectrum leakage 
problems and improve the detection accuracy of the harmonic 
parameters, but can not detect integer asked near the 
harmonic harmonic; cosine combination of high-end window- 
based dual-line[5,7,12] or line[13-14] interpolation FFT 
algorithm to estimate fundamental and the harmonic 
parameters, need to solve high-order equation[15-17], the 
computational complexity; continuous wavelet transform[18- 
19] can be realized between / harmonic detection, but the 
wavelet functions of different scales exist in the frequency 
domain interference, when the detection signal the harmonic 
components with frequencies close to, the detection method 
will fail; Prony method[20-21] is harmonic, harmonic 
analysis and modeling inter-effective way to accurately 
estimate the sinusoidal component of frequency, amplitude 
and phase angle , but the need to solve two equations and a 



polynomial of odd, high computational complexity and 
sensitivity to noise; there are other methods [22-24], or limited 
frequency resolution, or calculate volume, both in the specific 
application limitations. 

This paper presents a time-frequency filter, axle with the 
Hilbert transform in the frequency of negative, high-precision 
time domain convolution signal is detected between all the 
harmonics and the harmonic frequency, amplitude and phase. 
In this paper, a theoretical analysis and calculation formula is 
derived, the method to avoid the Fourier (FFT) domain 
spectral leakage, the entire sub-barrier effect and non-wave 
phenomenon. The simulation results show that: time- 
frequency convolution filter design flexible, easy to use, this 
algorithm can eliminate the harmonic interference and 
improve signal analysis precision, high accuracy for 
harmonic analysis. 

II. time-frequency filter design and analysis of 
Hilbert transform 

A. Time-frequency filter design 

Time-frequency filter: 

g(t, a> ) = sin(^— )e m >' (s(t) - s(t )) 

2 a (1) 

Where £ ^ is the step function, a ~ as ^ ne filter 
center parameters, the coefficient B is used to adjust the filter 
bandwidth (such as taking B = 1), coO center frequency. The 
frequency domain expression is: 
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Figure 1 Time-frequency filter of time / frequency characteristics 

Where 5a(f) = — p> & (0) = 1 figure 1 shows the trend of 
time-frequency filter characteristics, (a) trends in the time 
domain graph, (b) trends in the frequency domain; they 
change with the center frequency coO. By (2) and Figure 1 (b) 
shows G (co, coO) only in a narrow band centered coO 
significant amplitude, the other is almost zero. 

B.Hilbert transform 

HT is defined as the real signal [25] : 



f(t) = f(t)® 
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dr 
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Which ® for the sake of convolution, the Fourier 
transform F(co) , fit) the Fourier transform 

Pico) = -jSgn(co)F(co) (4) 
Sgn(co) the sign function. Signal f (t) through the 

Hilbert transform and its complex signal link, it can be 
analytic signal 

x(t) = f(t) + jf(t) (5) 
An important property of analytic signal is the positive 
frequency components retained, excluding the negative 
frequency part. That 

2F(ft>) (o > 

F{co) co = (6) 
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III Theoretical analysis and calculation 
formulas 

A.Analysis of Continuous 

If coO centered within the range of narrow-band frequency 
col>0 of the harmonic signal: 

f{t) = Acosicoj + (/)) 

Its frequency domain expression: 

Fico) = Anidico-co^e' 4, +S(a> + G> 1 )e- J '] (8) 

The analytic signal is the Fourier transform (from (6), (8), we 
have) 

Xico) = 2A7t8ico-co,)e 1 * 
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Y(a>) = G(a>,a> )X(a>) = H(co-a> )X(a>) 

= 2AnH{a> x - G> (t )e 1 ' 1 ' 8{a> - <j) x ) 
yit) = xit) ®git,co Q ) = AHico, - co Q )e j *e jm,t 



then: 
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The remainder mod[.] is divisible . 
B.Calculation of Discrete 

Set of discrete sampling frequency fs, the sampling 
1 

period DT= — ,Number of samples is N; 
fs 

take N 1 =[0.5N],N 2 =[0.94N],Computing discrete convolution: 

min{k-l,N] 

yik) = DT ^g(i)x(k-i) k = \,2,-,N (15) 

i'=max{l,i-iV} 

eik) = Argiyik))-Argiyik-\)) 

+ 2frmax{0, Sgni-Argiyik)) + Argiyik -1))} (16) 

k = N 1 ,---,N 2 ; Sgn is sign function 
The harmonic frequency f (Hz), amplitude A, the initial 
phase cp (°C) as, respectively: 

r AT, 

fs 



f = J\0(k) 

2niN 2 -N x +l)h, 

2ff_ 
fs 
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A = Y\y(k)\ (19) 

\H(co l -co Q )\(N 2 -N 1 +l) k % 

yr = (]£ {Arg( y(Jt)) + 2?r max (0, Sgn(-Arg(.y(jk)) ) - modl^t- 1 ^- ,2*]}) (20) 
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C. Experimental Evaluation 

Signal contains fundamental, DC, between 2 and 3 harmonic 
harmonic, and their parameters in Table 1, the expression: 

6 

f(t) = ^A i cos(a> i t + ti) 

'=0 (22) 
Its sampling frequency fs = 2kHz, the number of samples 
N = 5000. This method results are in Table 1 the right 
department. To the harmonic frequency, amplitude, initial 
phase of testing the value of the real value and are plotted in 
the same plot, the result is very accurate. 
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Figure 2 The harmonic characteristics of the true value of 
parameters compared with the measured values 

Table 1 harmonic of the frequency f = 350.7Hz specific icon 
near the detection algorithms 3, Figure (a), (b), (c) of the 
abscissa as the sample points. Figure (a) to (15) the magnitude, 
Figure (c) signal after filtering in frequency domain inverse 
Fourier transform (IFFT) of the amplitude, both in the same 
500 points; Figure (b) the type (16 ) transient frequency, 
Figure (d) the abscissa is the frequency (Hz), Hilbert 
transform signal was analytical signal in the frequency domain 
amplitude 
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Figure 3 the results of harmonic frequency f = 350.7Hz 

IV. CONCLUSIONS 

This new vision through the time-domain convolution can 
accurately detect the signal between all the harmonics and the 
harmonic frequency, amplitude and phase. In this paper, a 
time-frequency filter is designed with the Hilbert transform in 
the frequency of negative axle. The theoretical design analysis 
and calculation formula is derived, the method can be to avoid 
the Fourier (FFT) domain spectral leakage, barrier effect and 
non-integer secondary wave phenomenon. The simulation 
results show that: time-frequency convolution filter design is 
flexible, project implementation is easy and the algorithm is 
simple, response is quick. This algorithm can eliminate the 
harmonic interference and improve signal analysis precision, 
which can been applied in high accuracy for harmonic 
analysis. Applications of time-frequency domain filter 
selection bandwidth should be non-zero values within the used 
data segment, not to destroy its integrity 
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